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p ' Quantum Compiling Algorithms decompose (exactly, without approximations) an 

arbitrary 2^-^ unitary matrix acting on Nb qubits, into a sequence of elementary 



operations (SEO). There are many possible ways of decomposing a unitary matrix 



^ , into a SEO, and some of these decompositions have shorter length (are more efficient) 

than others. Finding an optimum (shortest) decomposition is a very hard task, and 
is not our intention here. A less ambitious, more doable task is to find methods 
^ ' for optimizing small segments of a SEO. Call these methods piecewise optimizations. 

■ Piecewise optimizations involve replacing a small quantum circuit by an equivalent 

one with fewer CNOTs. Two circuits are said to be equivalent if one of them multiplied 
by some external local operations equals the other. This equivalence relation between 
circuits has its own class functions, which we call circuit invariants. Dressed CNOTs 
are a simple yet very useful generalization of standard CNOTs. After discussing 
circuit invariants and dressed CNOTs, we give some methods for simplifying 2-qubit 
and 3-qubit circuits. We include with this paper software (written in Octave/Matlab) 
that checks many of the algorithms proposed in the paper. 
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1 Introduction 

Quantum Compiling Algorithms decompose (exactly, without approximations) an 
arbitrary 2^^ unitary matrix acting on Nb qubits, into a sequence of elementary op- 
erations (SEO). By elementary operations we mean operations that act on only a few 
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(usually 1 or 2) qubits (for example, all single-qubit rotations and CNOTs.) The most 
efficient quantum compiling algorithms to date are based on a recursive application 
of the Cosine-Sine Decomposition (CSD), a technique ffist proposed in Ref. IJ. An 
implementation of the algorithm of Ref. PP may be found in the computer program 
called Qubiter (patented, C++ source code publicly available at www.ar-tiste.com). 
Long after Ref.pj and Qubiter came out, many papers on quantum compiling via re- 
cursive CSD have appeared. These can be easily tracked down by making a keyword 
search in ArXiv or Google for something like ("Cosine-Sine" and "Decomposition" 
and "quantum") . 

We will call the number of CNOTs in a SEO its length. (Single-qubit rotations 
are not counted because these can be performed much faster than CNOTs.) Of course, 
there are many possible ways of decomposing a unitary matrix into a SEO, and some 
of these decompositions have shorter length (are more efficient) than others. The 
algorithm of Ref.jT] per se does not yield the shortest SEO. Finding an optimum 
(shortest) decomposition is a very hard task, and is not our intention here. A less 
ambitious, more doable task is to find methods for optimizing small segments of a 
SEO. Call these methods piecewise optimizations. The hope is that given any SEO, 
one can apply piecewise optimization methods to reduce the original SEO into an 
equivalent SEO whose length is much less, and might even be close to the shortest 
possible length. An analogy to our piecewise optimization strategy is the following. 
Think of a SEO as being like a path between two points in a manifold. If this path 
is initially unnecessary long, one might hope to make it a little less so by breaking 
it into pieces and optimizing the length of each piece. Breaking it into pieces again, 
and optimizing each piece again. And so on. 

Piecewise optimizations involve replacing a small quantum circuit by an equiv- 
alent one with fewer CNOTs. Two circuits are said to be equivalent if one of them 
multiplied by some external local operations equals the other. By external local oper- 
ations, we mean single-qubit rotations applied at the beginning or end of the circuit. 
This equivalence relation between circuits has its own class functions, which we call 
circuit invariants. Many excellent papers already exist on the use of such invariants in 
quantum computing. See, for example, Refs. [2], [S], and jH]. Such invariants are a 
crucial ingredient of this paper. (However, the paper does not assume that the reader 
possesses any prior knowledge about these invariants. The paper is self-contained in 
this regard.) 

Besides circuit invariants, another important ingredient of this paper is what 
we call dressed CNOTs (DC-NOTs). DC-NOTs are a simple yet very useful general- 
ization of standard CNOTs. To my knowledge, this paper is the ffist one to consider 
DC-CNOTs. DC-NOTs are convenient because they lump together a CNOT and 
some single-qubit rotations. Modulo external local operations, one can express any 
circuit solely in terms of a single type of circuit element (DC-CNOTs), rather than 
having to express it with two different types of circuit elements (CNOTs and single- 
qubit rotations). 
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After discussing circuit invariants and DC-NOTs, this paper gives some meth- 
ods for simphfying 2-qubit and 3-qubit circuits. 

Much is aheady known about simphfying 2-qubit circuits. Ref. jH] shows, via 
Cartan's KAK decomposition j7|, that a 2-qubit circuit with any number of CNOTs 
can always be reduced to a circuit with 3 CNOTs. Refs. [HI and [S] give necessary and 
sufficient conditions for when a 2-qubit circuit with 3 CNOTs reduces to fewer than 3 
CNOTs. In this paper, we spend some time re-proving these aheady known 2-qubit 
results using the new language of circuit invariants and DC-NOTs. This exercise 
yields new techniques and new geometrical insights that were lacking in previous 
proofs. 

In this paper, we also present some interesting new ways of simplifying 3-qubit 
circuits. Our results for 3-qubit circuits rely heavily on our results for 2-qubit circuits. 

We include with this paper software (written in Octave/Matlab) that checks 
many of the algorithms proposed in the paper. In the header of each section, and in 
the Table of Contents, each section name is followed by a list in square brackets of 
the names of the Octave m-files relevant to that section. Our software is not intended 
to be very efficient, or to be free of all conceivable loopholes. It is only intended to 
be a proof of principle of our algorithms. 

2 Notation 

[ global_declarations .m, global_def s .m, simul_real_svd.m, Gcimiiiajrep . m , 
sig.m, clieck_dciiots .m, f actor_SU2pow2_inatrix.m, f actor_SU2pow3jnatrix .m, 
test_f actor_su2pow.ni, get_normal_uiiit_vec .m, get_uiiit_vec .m ] 

In this section, we discuss notation, linguistic idiosyncrasies and abbreviations that 
will be used in subsequent sections. If any notation in this paper remains unclear 
to the reader after reading this section, he should consult Ref.|H|, a review article, 
written by the author of this paper, that uses the same notation as this paper. 

We will often use the symbol Nb = 0,1,2,... for number of bits, and Ns = 2^^ 
for the corresponding number of states. 

We will often abbreviate cos(a) and sin(a) by and Sa, respectively. We will 
often use a subscript of / to denote the final value of quantity that changes (e.g., a 
changes to df). When we say b = ±a, we mean b G {a, —a}. When we write X^^^, 
we mean, the quantity obtained by replacing a hj (3 everywhere in X. Likewise, 
by Xa^f^ we mean, the quantity obtained by swapping a and P everywhere in X. 
When we say ''A{ditto, A') is B{ditto, B'f we mean "A is B and A' is 5"'. LHS and 
RHS will stand for left-hand side and right-hand side. "RHON basis" will stand for 
" right-handed orthonormal basis" . 

Let Bool = {0, 1}. Let M denote the real numbers, C the complex numbers, Z 
all integers (positive and negative). For integers a and b, Tjaf> will denote all integers 
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from a to b, including a and b. If is anyone of the symbols >,>,<,<, and S is any 
set, define S^^ = {x E S : x Q 0} if the right hand side is defined. For example, 
are the positive integers. As usual, for any set S and r,p,q G will denote the 

set of r-tuples of S, and 5'*'^^, the set of p x q matrices with entries in S. 

As usual, U{Ns) will denote the Ns x Ns unitary matrices, and SU{Ns) the 

special (i.e., with determinant=l) elements of U{Ns). Given any A G U{Ns), we 

1 

define Ahy A = A/[det{A)] '^s , where we choose the principal branch of the function 

1 

(•)^s. We will refer to A as the "special counterpart" of A. (here the adjective 
"special" again means "with determinant=l"). 

will denote the set of all 3 dimensional real vectors, and = {a; G : 
\x\ = 1}. As is common in the Physics literature, a letter with an arrow (ditto, caret) 
over it, as in a (ditto, a) will denote an element of (ditto, W). a and a will be 
treated as column vectors when they appear in matrix expressions. 

Let Sj G M.^ for j G Zi^^- We will use the following non-standard notation for 
r-fold cross products: 

[010203 . . . Or) = (• • • ((oi X 02) X 03) • • • X Or) . (1) 

For example, [01020304) = ((oi x 02) x 03) x 04. Of course, an (r-|-2)-fold cross-product 
can be reduced to an r-fold cross-product using the well known "BAG minus GAB" 
identities: for o, 6, c G IR , ax{bxc) — b{a-c)—c{a-b) and {axb)xc— b{a-c) — a{b-c). 
For example, if o, b are perpendicular unit vectors, then [abb) = —a. 

Suppose o, 6 G M^. angle{a, b) will denote the angle between a and b, defined 
up to 271. We will say a is parallel to b and write o || 6 iff o x 6 = 0; i.e., iff o = ±6, 
or o = 0, or 6 = 0. We will say a is perpendicular to b and write o _L 6 iff o • 6 = 0. 
For 6 7^ 0, define Oyg, the part of a along 6, by 

— * — * 

For 6 7^ 0, define a^j, the part of a across 6, by 

^ ^ ^ {d-b)b -[abb) 
a.^ — a — Onr —a — = — - — . (3) 

\\b |^|2 ^ ' 

For any square matrix A, A^ will denote its transpose, A* , its complex con- 
jugate, and A^ = A*'^, its Hermitian conjugate. Sij will denote the Kronecker delta 
function. (It equals one if i — j and zero otherwise.) 

Let l2,crx,crY, Oz be the 2d identity matrix and Pauli matrices. Sometimes, we 
set (Xi, X2, X3) = (X, F, Z) and denote the Pauli matrices by ox^, (yx-^-, ^Xz- Suppose 
W G {X, F, Z}. Define the number operators: nw = ^-f^ and nyv = Note 
that (— 1)"'W' — Usually, nz is denoted merely by n and nz by n. If Wj G 
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{1,X,Y,Z} for j e Zi^Nb, let awi,W2,...,WMg = <^Wx ® (yw2 ® ■■■f^M/jv^> where any 
incidence of ui on the RHS is replaced by 12- For example, (Txyi = cx ^ cry ^ h- 

H — ^ is the one-bit Hadamard matrix and H^^'^ is its A^B-fold 

tensor product. H satisfies H'^ — 1, HaxH — az, HazH — ax and HayH — —ay- 
Suppose ao G M and a G M^. Wc will abbreviate c? ■ a by ag. The standard 
terminology is to call + ^ quaternion, and to call a^ a vector quaternion 
(divided by i). To shorten this terminology, we will refer to ag as a Paulion, and call 
a its defining vector. If |a| = 1, we will call as a unit Paulion. One can reduce a 
product of two Paulions by using the identity crgcrg — a-b + For a e R^, define 

number operators na = and fia = Note that (-l)"" = a^. If Wj e W or 

Wj = 1 for j e Zi^Nb, let awi,W2,...,WNg = (^Wi ® o^Wi ® • • • c^tyivs ' 

Suppose Ai is the set of all matrices M G C'*^'* that can be expressed in the 
form M = Y^f, a^^ g^, where a^, G for all k. Suppose C is the set of all matrices 

L G M^^^ that can be expressed in the form L = Ofc^fc , where dk,bk G for 
all k. For every M G let r(M) or represent the 3x3 matrix with entries 
jtT:{aXi,XjM), where i,j G Zi^^. (The symbol F was chosen to evoke the mental 
picture of a column vector times a row vector; such is the output of the function 
r(-)). For every L E C, define r~^(L) = j o'x,,XjLij. It's easy to check that 
rr~^ = r~^r = l so the map r : — * £ is 1-1 onto. Let lin{Ai) be the set of 
linear combinations over C of elements of Ai, and lin{C) of C The map F can be 
extended to F : C + lin{M) ^ C + lin{C), T{X + ^. aiMi) = A + a^Mf . T js 
also a 1-1 onto map. Henceforth, we will use F to refer to both F and its extension F. 
Given a matrix A G C + lin{Ai), we will call its Gamma representation. Often, 
in contexts where this will not lead to confusion, we will drop the F superscript and 
denote simply by A. 

The next theorem, although almost trivial, will be used frequently in this 

paper. 

Theorem 1 The map / : x M'^ — SU{2), f{a,b) = aaCf, is well defined and 
onto. In other words: (well-defined) If a,b G M^, then f{a,b) G SU{2). (onto) If 
U G SU (2), then there exist a, 6 G such that U = f{d, b). 

proof: 

(well defined) Given a, 6 G M^, one can always find an angle 9 such that d-b = cq 
and la X 61 = sq. Let w = It follows that a^at = d-b + ia-^y^ = e'^"^ G SU(2). 

(onto) Given U = e*^'^*, where w G M-^ and G M, one can always find a 
(non-unique) pair of unit vectors d and b in the plane perpendicular to w, such that 
9 — angle{d, b), and dxb points in the w direction. Hence, d-b — cg and dxb — sgw. 
It follows that (7a cr^ — d -1) + icr-^^ = g2^o-ta 
QED 
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One has: 



afaaCTf = Crf(0"a||^ + C^d^f)(^r = (^d^^f. " (^a^f. = 0"a||,,-axf = (^df ■ (4) 

A geometrical interpretation of this identity is shown in Fig^. The similarity trans- 
formation af{-)af takes the Paulion aa to ua^, where a/ is the reflection of a on f. 

Suppose d,b E M.^, and we want to find U E SU (2) such that erg = WaaU. 
Such a U can be constructed as a product of two Paulions (See FiglDo). Indeed, let 
9 = angle{a,b) and p = jr^. Let f be the vector that bisects the angle between a 

and b, and is oriented so that axf points along p. Note that b can be obtained by 
reflecting a on the bisector f. Hence 

aaCTf = e'^^'P , O-g = afaaCTr . (5) 

Combining these two results yields 

cxg = (a^(Ta)aa(fTa(Tf) = e'^^'^Vae^^'^^ . (6) 



A 




(a) (b) 

Figure 1: (a) If arCraCf- = o^dj, then a/ is obtained by reflecting a on f.(b) Suppose 6 
is the result of rotating a by an angle 9. Then 6 can be obtained by reflecting a on 
the bisector f of the angle between a and b. 



3 Invariants for Quantum Circuits 

In this section, we will discuss circuit invariants; i.e., functions that map all equivalent 
circuits to the same value. By equivalent circuits we mean circuits that are equal, 
modulo external local operations. 

Suppose A and B are elements of U{Ns) ( i.e., they are A^'^-qubit gates). We 
will say A and B are equivalent under local operations on the right hand side 
(LO-RHS), and write A ~/j B, iff there exist Uj E U{2) for j E Zq ^b-i such that 

B = A{Unb-i 0...0U2 0Ui®Uo) . (7) 
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is clearly an equivalence relation as it is symmetric, reflexive and transitive. 

Henceforth, we will say that a function x with domain U{Ns) is a LO-RHS 
invariant if for any A,B E U{Ns), A B implies that x(^) = for some 

C G M (C may depend on A, B). 

A frequent goal is to find a complete set of scalar invariant functions; that 
is, a set of functions Xj '■ U{Ns) — > M such that for any A,Be U{Ns), A ~r B if 
and only if Xj{^) = Xj{B) for all j. An extensive literature already exist on such 
invariants. They were first studied by Group Theorists, and, in more recent times, 
they have been used by Quantum Computerists [2], [S], 

One can define an analogous equivalence relation for local operations on 
the left hand side (LO-LHS), and an equivalence relation r^LR for local operations on 
both sides (L0-2S). Of course, the equivalence classes (e-classes) of are a disjoint 
partition of U{Ns). Ditto for the e-classes of ~i and ^lr- It's also clear that any 
e-class for ~^ is contained in an e-class for ^lr, and that some e-classes of ^lr 
contain more than one e-class of (In fact, the e-classes of contained within a 
single e-class of ~Li?, can be labeled by the elements of f/(2)®^^). 

Note that for any a G M^, 

ayajo-y = -0-5 . (8) 

Hence, for ^ G M and a G M^, 

ay[e'('^+^^Y^Y = e*^'-"'^^ • (9) 

Thus, when U G SU{2) (but not when U G f/(2)), ayU^ay = U'^ = W. 

For any A G U{Ns), define a quadratic (second order in A) invariant 

= Aay^^M^ay®^^ . (10) 
For example, for A G f/(4), A^^) = AayyA'^aYY- 

Theorem 2 

(a) For A,B e SU{4), A^rB if and only if A^^^ = (-l)"^^^) for some n E Z. 
(h) For A,B e U{4), Ar^RB if and only if A^^) = e^'^E^^) for some C G M. 

proof: 

(a) Assume A,B E SU{4). A can always be represented in the form 

A = r^^^ exp{iajk(Tx,Xk) exp(ia'.(Tx,i) exp(iafccrixj , (H) 
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where n{A) G Z and ajk, a'j, ak G R. (Note that det(z/4) = 1 so det(A) = 1.) We 
are using Einstein's imphcit summation convention, and j, k range over {1, 2, 3}. 
By Eqs.dHI) and (HD), 

aYvA^aYY = exp(-mfc(JixJ exp{-ia'jax,i) exp{iajkaxjxj ■ (12) 

Thus 

A(') = (-l)r^(^) exp{t2a,kCTx,x,) ■ (13) 
Likewise, B can be represented in the form 

B = r(^) exp{ibjk(rx,xj exp{ib'jax,i) exp{ibk(TixJ , (14) 
where n{B) G Z and 6^-, bk G M. Then 

5(2) ^ (_i)n{B) exp(z26,fcax,xj ■ (15) 

(^) Suppose A ~R B. Looking at Eqs.©, (fTTj) and (HH), we see that for every 
j, k, there exists an integer njk such that ajfc = bjk + nrijk. Therefore, 

exp{i2ajk(rx,xj = exp{i2bjkax^xj ■ (16) 

Therefore, looking at Eqs. pHj) and (fTH|) . we see that there exists an integer n 
such that = (-1)"5(2). 

(<^) Suppose = (-1)'^5(2). Then, looking at Eqs.jni) and we see 

that for every j, k, there exists an integer njk such that 2ajk = 2bjk + vrrzjfc. 
Therefore, 

exp(ia^fcax,xj = exp{tbjk(rx,x,) Ylii^^x.x.P' ■ (17) 

Therefore, from Eqs.©, (fTT|) and (HH), we see that A ~r 5. 

(b) Assume A,B G f/(4). Eqs. lfTTj) and (fT^ still apply except that we must replace 
in them by e^'^(^) and (-1)"(^) by e*2C(^) for some ((A) G M. Eqs.dH)) and 
fll5|) still apply except that we must replace in them by e*''*-^-' and (— 1)"-(-^) 
by e*2C(^) for some ({B) G R. 

(^) Suppose A r^ji B. Eq. lfT^ still applies so there exists C ^ R such that 
(^) Suppose = e*^S(2). Eq.(|I7D still applies so A ~r 5. 
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QED 

By virtue of Theorem |21 the absolute value of the entries of the matrix A^'^^ 
are a complete set of LO-RHS scalar invariants for Nb = 2. Theorem |2fa) reflects 
the fact that when A,Be SU{4), since A and B must both have unit determinant, 
the only local operations connecting A and B are either elements of SU{2) or i or 
products of these. Applying an SU{2) gate to the RHS of A does not change A^'^\ 
whereas applying i changes A^"^^ to its negative. 

Now suppose Nb = 3. One can represent any A G SU{8) as 

A = e'^""^^^ exp{iajkr(TxjX,Xr) 

eyip{ia'-,,aiXjxJ exp(m^;,(Tx,ixJ exp^iajkaxjXki) 
exp{ia'-ax^u) exp{ia'jaixji) exp{iajauxj) ■ (18) 

When the continuous parameters of A are small, 

^(2) ^ g*in{A)[^ ^ 2t{a%aix,x, + a]^ax,ix, + %fc^x,x,i)] • (19) 

This A^'^'^ is independent of the ajkr parameters. So, for A, B E SU{8), A^'^^ = ii?*^^) 
or A^'^^ = iiB^"^^ is a necessary but not a sufficient condition for A B. More 
invariants than just A*^^) are needed for Nb > 2. 

Higher order invariants can be generated as follows. We will represent them 
diagrammatically using the symbols defined in Fig|21 FigOl shows second and fourth 
order invariants under LO-RHS for a circuit with 3 bits. The same idea can be used to 
generate invariants of order equal to any even number, for any number of qubits. FigH 
explains why the circuits portrayed in FigI21 are invariant under LO-RHS. Roughly 
speaking, if we apply aU E SU{2) to the RHS of A G SU{8), then, in the diagram 
of a fourth order invariant, a copy of U must be inserted next to each of the 4 copies 
of A. And these 4 copies of U annihilate each other. This paper will only use the 
second order invariant A^^\ We will not even use Group Theory in this paper. For 
information on the group theoretic underpinnings of quantum circuit invariants, see, 
for example, Ref.[2]. 

Figure 2: Key to symbols used in Figs 01 and El A G SU{Ns). 
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Figure 3: Second and fourth order invariants under LO-RHS for a circuit with 3 bits. 
A G S[/(8). 



4 Dressed CNOTs 

[ drll.m, drllO.m, drOll.m, drlOl.m, ] 

In this section, we define dressed CNOTs, a simple yet powerful generalization of the 
standard CNOT. We also discuss some simple properties of dressed CNOTs that will 
be used in subsequent sections. 

The controlled NOT (CNOT) with control bit 1 and target bit 0, is defined 

by 

= _ ^^(o)"(i) . (20) 

Now suppose U and V are arbitrary elements of SU{2). Define a and a! by UaxU'^ = 
a a and VozV^ = cr^'. Then a dressed CNOT (DC-NOT) connecting bits and 
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Figure 4: Why the diagrams shown in Figl^lare invariant under LO-RHS. A G SU{8) 
and U e SU{2). 



1, is defined by 




— V 



(_l)"a(o)na'{i) = (Ta(0)"a'(^) = • (21) 



We will refer to the vectors d' and a as the defining vectors of the DC-NOT. 

Sometimes in this paper, we will draw a circuit containing one or more DC- 
NOTs whose oval nodes are empty. By this we will mean that the omitted defining 
vectors are arbitrary and their precise value is unimportant in that context. 

Consider the wire corresponding to bit /i in a quantum circuit. Within the 
bit-/i wire, consider two adjacent oval nodes belonging to two different DC-NOTs: 
— (a) — (bj — . If 6 II a, we will say there is a breach at that position in the bit-// 

wire. If b -L d, we will say there is a foil at that position in the bit-/i wire. 



Theorem 3 




-(1 + 0"l,a + CTa'A - (Ta',a) 



proof: 



(22) 



a'(l)"-(°) = aa'(l)ria(0) + rra(0) 



(23) 



QED 
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proof: 

aa(0)2"a'(i) = 1. 

QED 

Theorem 5 




(25) 



proof: 

i-aaiOT"^'^ = (-l)"^'('Va(0)"'^'(i) = aa'(l)aa(0)"^'(i) . (26) 

QED 

In subsequent sections, we will often need to calculate the effect of a similarity 
transformation produced by pre and post multiplying an operator by the same DC- 
NOT. The next theorem will be useful for performing such calculations. 

Theorem 6 



proof: 

Clearly, 




On the other hand. 




QED 

5 Wake Identities 

In this section we prove what we call a "wake identity". We call it thus because in it, 
one DC-NOT is pushed through another, producing a third DC-NOT as its "wake" . 
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Theorem 7 Suppose a! A-V . 



proof: 




QED 



6 Swapper Identities 

[ swap_t3.m, test_swap_t3 .m ] 

In this section, we discuss certain DC-NOT identities associated with the qubit Ex- 
change Operator (a.k.a. Swap Operator or Swapper). 

We will represent the Swapper by a double arrow connecting the two qubits 
being swapped. By definition, the Swapper satisfies 



U 



-Ttr- 



-V- 



^ U 



(32) 



14 



for any U G U{2). As is well known (for a proof, see, for example, Ref.|H]), the 
Swapper can be expressed clS cL product of 3 CNOTs: 



nz: ■ 



(33) 



The next theorem shows that the Swapper can also be expressed as a product of 3 
DC-NOTs. 

Theorem 8 Suppose a -Lb, U E SU{2), WaJJ = Oo; , and U'^aylJ = ay. 






\a) \h 

161 la 




U 



(34a) 



(34b) 



proof: 



Then 



Since a _L 6, there exists V G SU{2) such that V^axV = Oa and V'^azV = erg. 



-Ttr- 



-A 



-V- 



— — Ft 



V — 



V 




■ (35) 



This proves Eg. ()34a|l . Eg. ()34bjl follows from 

4 



f/ — 



-V- 



u 



(36) 



QED 

We will refer to the next identity, Eg.(jnZj), as the 2/3-Swapper identity, because 
its LHS contains 2/3 of a Swapper. 

Theorem 9 For any a G M, 






Jixy 


) — ii) — 


U 








) — v£) — 


U' 






u 


— w' — ( 








U' 


— ^ — ( 


Pzy^ 





(37a) 



(37b) 
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Figure 5: Orientation of vectors o' , and i6' ,• Note that Ha^i H = a^' ■ The same 
picture, but omitting all primes, describes q^y and 



where (see Fig\^ 



Pzy CaZ + Sai) , P^y iPzy)a 



{3t 



Qxy CqiX ^aV ) Qxy i^Qxy^a—fa' 5 ("^Q) 

(p vector has a positive sign in front of s^, q vector has a negative one) and 

U = e'^'^^e-'^"^ , U' = (f/)aw . (40) 

Note that the left-hand sides of Eqs. f37a}) and l{37b\ ) are independent of the two angles 
a and a! ; only their right-hand sides depend on these angles. 



proof: 



From the expression of Swapper as a product of 3 CNOTs, we get 




E 



H — 



(41) 



From FigEl it follows that 



(42) 
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Thus 
















2 g * 2 


0-x 



(43a) 



(43b) 



(43c) 



(43d) 



QED 



The next theorem follows immediately from the previous one, by a change of 



basis. 

Theorem 10 Suppose a e a ± 6, and a! A-V . Then 





(44a) 



(44b) 



where 



O/ = CqO + Sa.[(ih) , a'j: — Ca'd! + Sa'[a'b') , 
bf ^Cab- Sa[ab) , b'f = c^'b' - Sa'[a'b') , 



and 



■ a _ -a 



(45) 
(46) 

(47) 



17 



proof: 

Just change basis in the space where bit (ditto, bit 1) hves so that (x, z) 
is replaced by {b, [ab),a) (ditto, (6', [d'b'),a')). 
QED 

We will refer to the next identity, Eq. ipHj) . as the 1/3 Swapper identity. 



Theorem 11 





(4^ 



f/'t- 



where all variables are defined as in Theorem\^ 
proof: 



From the Hermitian conjugate of Eq. ()37a|) . one gets 

(Qxy) 





(49) 



Let LHS and RHS stand for the left and right hand sides of Eq.(g8j). Pre- multiplying 

both sides of the last equation by T yields 

iz) 



LHS 




RHS . 



(50) 



QED 



7 DC-NOT Similarity Transformation Identities 

[ sim_trans_t4 .m, test_siin_trcLns_t4 .m ] 

In this section, we present some identities which contain a similarity transformation 
produced by pre and post multiplying an operator by the same DC-NOT. 

We will refer to the next theorem as the DC-NOT similarity transformation 
identity. 

Theorem 12 For any a, A G M, 




Cgax + Sgaz \ —{x)— 



(51) 
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where qxy = c\x — s\y. Note that the LHS of Eg. \51]) equals its RHS evaluated at 
A = 0. 



proof: 



Since 



it follows that 



[qxyz) = -{cxy + SxX) 



(52) 



^qxytixy + ^[qxyZ)[qxyz) <^ CxX- S xV ,CxX- S xV + CxV+S xX ,Cxy+S xX (53a) 

(53b) 



'^xx T O'i/y ■ 



Let LHS and RHS denote the left-hand side and right-hand side, respectively, of 
Eq.(jniD. Then, using Ea. (|77j) . 



(Qxy ) \M^y ) ylxy ) yl^y J 

RHS = (SaO-q^y,! + CaCr^^i) {(^aC^lAxy + Sa(Tl,z) Jl54a) 

[Qxy J [Qxy J [Qxy J 

= SaCaicTg^yg^y + cr[g^yz)[qxyz)) + c^a ^1 + s^a u 

= SaCaicTrtx + CTyy) + cl^CTzl + s1^<^iz 



QED 



LHS . 



It is convenient to define, for any ^ G 




(55) 



(The p vectors have a positive sign in front of the sine function whereas the q vectors 
have a negative one). 

The next theorem follows from the DC-NOT similarity transformation identity. 



Theorem 13 For any 0, A G M, 

{Pzx 
U 

where (see Fig[^ 






U 


— fa'A — 


U' 



(56) 



(57) 
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and 



where 





^out 



Figure 6: Variables used in Theorem IT^ 



U = {CaCFx + SaCrz)(Ca(T^A + Sa(Tz) , U' = (U) 



(5^ 



TT 



2a = - - , 2p = n-2a . 



(59) 



proof: 



From Figini it follows that 



and 



(60) 



(61) 



Let LHS and RHS denote the left-hand side and right-hand side, respectively, of 
Eq.dSni)- Then 



LHS 



e 


iaay 




e 


ifScry 




(62) 
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and 



RHS 



Therefore, Eq.(|3Bj) is equivalent to the assertion that 






Now pre-multiply each side of the last equation by axx 



(63) 



(64) 




Ca<7X + SaCTz 



C/3crx + spaz 





xyl 



(65) 



The preceding equation follows from Theorem IT^ and the fact that a + (3 = tt /2. 
QED 

The next theorem is a simple variation of the previous one. (The left-hand 
sides of Eqs.(jHni) (jMj) differ only in that one circuit has two g's in the bit-1 wire 
whereas the other circuit has two p's.) 



Theorem 14 





(66) 



where all variables are defined as in TheoremM'A 
proof: 

Let LHS^^ represent the left-hand side of Eq. (|3^ . and LHS^^ the left-hand 
side of Eq. . Then 



— 0"aA ,crx — 



LH 



oz 



oz 
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The right-hand sides of Eqs. ()56p and ()66p must be related in the same way as their 

left-hand sides. 

QED 



8 LO-RHS Invariant for Circuits with 
Two Qubits, and Multiple DC-NOTs 

In previous sections we defined the LO-RHS invariant A^"^^ for any A G U{Ns). We 
also defined DC-NOTs and discussed some of their properties. In this section, we 
combine these two concepts: we calculate A^'^^ when A is a product of one or more 
DC-NOTs acting on the same two qubits. 

Henceforth, we will denote the product of r DC-NOTs (all acting on the same 
two qubits) by the symbol Qr followed by a list (enclosed in parenthesis) of its ar- 
guments. Sometimes, if this doesn't lead to confusion, its hst of arguments will be 
omitted. Thus, 



02 CLl 

d'2 a'l 




The determinant of Qr equals either plus or minus one. Indeed, 




(68) 




det 



det 



h 
ax 



Since det (AS) = det (A) det(S), it follows that for r = Z 

det(a) = (-1)' . 
It is convenient to define a matrix Qr by 



>0 



(69) 



(70) 



(71) 



Henceforth, we will refer to Qr as the special counterpart of Qr- (Here the adjective 
"special" means "having unit determinant"). Qr G f/(4) and det{Qr) = 1, so Qr G 
SU{A). 

Since aycrjay = o"_a. 





(72a) 
(72b) 



22 



For r G Z^°, Qr obeys the following recursion relation: 




Note that the LO-RHS invariants of Qr and of its special counterpart are 
related by 

Gi'^ = €Qf^ . (74) 

The remainder of Section |H1 consists of 4 subsections which give explicit formu- 
las for Qi- for r from 1 to 4. These 4 subsections are very useful, but make for dry 
reading when considered in isolation; they only come alive and prove their mettle as 
we start using them in subsequent sections. Thus, the reader is advised not to spend 
too much time on them during his first reading of this paper. He should skim the 4 
subsections, and then come back to them as the need arises. 

8.1 Invariant for Circuits with 1 DC-NOT 

[ ckt_invarl23.m ] 

This part of our program is dedicated to the letters Qf^ . 
Theorem 15 



proof: 




QED 

8.2 Invariant for Circuits with 2 DC-NOTs 

[ ckt_invarl23.m, diag_ckt_invar2 .m, diag_ckt_invar2_aux .m, 
test_diag_invar2.m ] 

(2) 

This part of our program is dedicated to the letters ■ 
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Theorem 16 



where 




Xor + i)\2i + Mr + i^' 



2i 



(77a) 



(77b) 



A2r = (o- ■ i){ci ■ h') 



(78) 



Asi = , 



(79) 



A2r — ^[a'yy)\ahh) ' 



(80) 



A2i = a ■ &(T^,xS',S + « ■ ^'^S',axS • (81) 

proof: 

fn\ 

An explicit expression for g\'> was given in Section 18.11 Eq. ()27p shows how to 
calculate the effect of DC-NOT similarity transformations. Using these two results, 
one gets 



QED 

Theorem 17 



\—cra'h\ 




{a! ■h'){a-h) 



-a 



a'b'b'),[dbb) 



[A2r,A2i] = 



(A^,)^A^, = , 




(-1) 



(82a) 

(82b) 
(82c) 

(82d) 



^3) 
^5) 
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This follows easily from Eqs. ()80|) and ()81|) . 



proof: 
QED 

It is convenient to parameterize the expression for Q'^^ given by Theorem 
using as few parameters as possible. 



A A 



/ = 




[ah) 



AAA 

[abb) 



lab'} 




'^1 '^z 



-AAA 

/ a b' b') 



(2) 

Figure 7: Principal parameters of ^2 ■ 



Theorem 18 Q2 can be parameterized with 2 real numbers a, a' , and 2 RHON bases 

{} j)j =1^2, z md {f'j)j =1^2, ■i- Call these the principal parameters of Q^^ (seeFig^. More 
explicitly, 

= \2r + A2r + l^2i , (86) 

where 

^2r = Ca'Ca , (87) 

A2r = -{Sa'Sa)fVl , (88) 



A. 



2i 





/r 


f7 









fi 















(89a) 
(89b) 





Co' Sa 



[h h] ■ (89c) 
(Eqs J89a\) . Ii89bp . and ^ESB) are 3 different styles of representing the same thing.) 
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proof: 

Define a' e [0, tt) to be the angle between a! and h' . Thus 



Cq,/ = a' • 6' , Sq' = |a' x h'\ . (90) 

If Sq/ 7^ 0, set 

If Sq,' = 0, choose (/j)j=i,2,3 to be any RHON basis with /( = h' . 

Use the previous paragraph with all primes removed to define a and (/j)j=i,2,3- 

QED 

Suppose we are given a matrix which is known to be the LO-RHS invariant 
Q^^ of a quantum circuit with 2-qubits and 2 DC-NOTs. Furthermore, we are asked 
to extract from this matrix values (non-unique ones) for a, 6, a' and h' . Next we 
will give an algorithm for accomplishing this task. We will call it our "Algorithm 
for Diagonalizing Q^^'' . The algorithm first expresses ^2^'* in term of its principal 
parameters. Then it solves for a, 6, a' and h' in terms of these parameters. 

Algorithm for Diagonalizing Q2^: 

1. Set \2r = -MQ?^). Set A = - As., A2. = (A + At)/2 and A2, = (A - 

A^) / (2i) . Hence, Qf^ — + ^2r + "1^211 where As. is a real scalar, and As. , Asj 
are traceless Hermitian matrices. 

2. Calculate Cq/Cq,, Sa'Sa, /2 and /s from As. and h.2r- (If ^2r — 0, then take 
Sa'Sa — 0, and choose /s and /s to be any 3d unit vectors.) 

3. Choose any RHON basis (/ij)j=i,2,3 such that = A, and any RHON basis 
{h'j)j=i,2,3 such that h'2 = /a. 

4. Find a Singular Value Decomposition (SVD) of the matrix 



M 



(92) 



h'^A2ihi h'^A2ih3 

In other words, find 2-dimensional orthogonal matrices U, V and a non-negative 
2-dimensional diagonal matrix D such that 



M = UDV^ . (93) 
Now calculate s^/Ca, Ca'Sa, /a, /i, /a, h from 



^Oi'^Q. 



= ^ , (94) 
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[flIi] = [KK]U, (95) 

and 

[fij3] = [hi,h]V . (96) 

5. By expressing U on the RHS of Eq. ()95p in component form, it is easy to verify 
that 

X /{ = det{U)h', X h[ . (97) 
/ig X h[ ■ h2 = +1 and = /I'a so 

/^x/(-/^ = det(f/). (98) 

det{U) will always equal either +1 or —1. If det(?7) = —1, replace /g — > — /g 
and Sa'Ca — — Sa/Ca. Thcsc replacements make (/(, /g, /g) a right handed basis. 

If det{V) = —1, an analogous procedure can be used to convert (/i, /2, /s) into 
a right-handed basis. 

6. At this point, we know the four quantities c^'Cq, Sa'Ca, Ca'Sa, and Sa'Sa- Calcu- 
late a' ± a from 

COs(q;' ± a) = Ca'Ca T Sa'Sa , (99a) 

and 

sin(Q;' ± a) = Sa'Ca ± CQ,/Sa . (99b) 
Calculate (a' a) from a' ±a. 



7. Calculate a, b, a', b' from: 



^ = A ^ . , ( ^' ^ " 

a = Ca/i - Saf2 ' 1 a' = Ca//{ - Sa'f2 

8.3 Invariant for Circuits with 3 DC-NOTs 

[ ckt_invarl23.m, ckt_invar3 .m, diag_ckt_invar3 .m, test_diag_invar3.iii ] 

(2) 

This part of our program is dedicated to the letters . 



(100) 
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Theorem 19 




where 



Asr = < 



Asr = [ab'b') ■ c [abb) ■ c , 
A3, = -(a ■ ■ c)V - (a ■ b'){b' ■ c')V , 
-(a' ■ b'){a ■ b)ac',c 

+ (a ■ ■ c)a^,,y,,) , + {a' ■ b')ib' ■ c')a,,^^,^,^ 
+ (a'-6')Va[j,£,),£+(a-6)VV,,[g,) 

~'^[a'yb'c'c'),[abbcc) 

where V = a x b ■ c and V = a' x b' ■ c'. 
proof: 



A 



3i 




(101a) 
(101b) 

(102) 
(103) 

(104) 
(105) 



(106a) 
(106b) 



An explicit expression for q!^^ was given in Section lST^ Eq. (j27p shows how to calculate 

the effect of DC-NOT similarity transformations. 

QED 

Theorem 20 Suppose 



C 




n 




(107) 



For any C, it is possible to find an TZ such that C ^rTZ, and such that (ajaj xbj -cj 
0, and (h)b'f ± span{c'j , a'jr) . 
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proof: 

As pointed out in the introduction, Ref. 6J shows how to express any 2-qubit 
unitary operation as a circuit with just 3-CNOTs. It is easy to check that conditions 
(a) and (b) are satisfied by the 3-CNOT circuit given in Ref.|n]. Hence, this theorem 
has already been proven in Ref. [H], although Ref. [HI does not explicitly point out this 
property of their 3-CNOT circuit. The "Algorithm for Diagonalizing ^3 , that is 
presented later in this section, also constitutes a (constructive) proof of this theorem. 
QED 

For A,B& RP^'J, define the following two commutators: 

[A, B]l = A^B - B'^A , (108a) 

[A, B]r = AB^ - BA^ . (108b) 

(Here, the letters L and R stand for left and right. They indicate on which matrix 
the transpose symbol acts, either the left or the right matrix in the matrix product.) 
Ref. [7] presents a proof (due to Eckart and Young) of the following Theorem. A,B& 
j^pxg j^g^yg slmultaueous Singular Value Decomposition (SVD) if and only if [A, B]l 
and [A,^]/^ are both zero. By a simultaneous SVD we mean orthogonal matrices 
U, V and real diagonal matrices Da, Db such that 

A = UDaV^ , B = UDbV^ . (109) 

When considering the SVD of a single matrix A, one usually insists in making the 
entries of Da non-negative, and calling them the singular values of A. In the case 
of a simultaneous SVD, one can't always make both diagonal matrices non-negative, 
but one can certainly make one of them so. 

Of course, the previous paragraph applies almost intact if A and B are elements 
of C^^"? instead of M^^'^. For A, B complex, one must replace the T (transpose) symbol 
by the f (Hermitian conjugate) symbol in Eqs. ()108p and ()109p . Also, the matrices 
[/, V in Eq. ()109p must be unitary instead of orthogonal. 

Note that when A and B are Hermitian, the condition that [A, B]l and [A, B]ji 
both vanish becomes simply the condition that A and B commute. The Eckart, Young 
theorem then becomes a theorem very familiar to practitioners of Quantum Mechan- 
ics: two Hermitian operators can be simultaneously diagonalized iff they commute. 

Theorem 21 

[A3r,Ad = 0, (110) 

[aLA3Jl = o, (111) 

[kl.,Kl]R = Q. (112) 
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proof: 

Let 



A = ^f -tr(^f), (113) 

so 

A + At ^ A- At ^ ^ 

As. = ^ , A. = ^— . (114) 

Thus, 

[A3., A3.] = ^[A + At, A - At] = i[At, A] = ^ ] = , (115) 

where the last commutator is zero because ^3 is unitary. 
Note that for any a, a', 6, 6' G M^, 

= m^-b)a,,,y^,-ia'.b')a,^,J. (116b) 
From Theorem ^1 we know that A^j. and A^^ can be expressed in the form 

Asr = "^ajaa'^^a, , A3i = ^ Pk(Tyj^ , (117) 

j k 

for some Oj, /3j G M and d^, 6^, 6'^ G M^. Therefore, 

= [A3.,A3i] (118a) 
= Xl^J^^^N'^^'^^'^sJ (118b) 

= z2 ^ aj(3k[{aj ■ bk)a^,^y^^^ - (a^ ■ b'^)a^^^^^f^J . (118c) 

This imphes that 

(yjPk{aj ■ bk)aj xb'i. = , ^ ajPkia'j ■ b'^)aj x bk = . (119) 
Now note that 

[A^r.A^R = (^a,a;aJ)(5^/?fcML^)-(5^/?fc6l6D(5Z«.«.«r) (120a) 

= Y,a,Pk{aJbk)[a'/b'^ -b'^a^^] (120b) 

= , (120c) 
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where the last expression vanishes due to Eq. ()119p . An analogous argument shows 

that [Ag^jAgJi also vanishes. 

QED 

It is convenient to parameterize the expression for Q^^ given by Theorem I19| 
using as few parameters as possible. 

A AA 



-[ b c c ) 




Figure 8: Principal parameters of 



(2) 

Theorem 22 ^3 can be parameterized with 3 real numbers /5,/?i,/?2; O'lT'd 2 RHON 
bases 2,3 O'^d (^j)j=i,2,3- Call these the principal parameters ofQ^^ (see Fig\^. 

More explicitly, 



where 



where 



= A3, + ^Aa. + A3,, + zAg, , (121) 

As. = -Xo , (122) 

Asi = -Yo , (123) 
3 



A3r. = (124) 

3 

A3. = (125) 



Xo = cp^sp^sp^ , (126) 
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(127) 



1 2 3 

2 3 1 



(l^i)j=l,2,3 = (S/3C;3iS/32, SpSp^\CfS,\,CpCfS,CfS,) , (128) 
(/^j)i=l,2,3 = (-C/jS^Jc^J, -C^Cfs.Sfs,, Sfs^S^^Sp^) , (129) 

where ^ G {+1, —1} one? 7r() is the permutation 
proof: 

We will assume from the onset of this proof that (a)d x b ■ c = 0, and {h)b' _L 
span{cf , a'). This can be assumed without loss of generality because of Theorem 1^ 
Let 

^ = sign([a6) ■ [be)) , ^2 = sign(6 ■ c) . (130) 

Without loss of generahty, we will assume that —^^2 = +1- If ~^^2 is initially 
negative, we can make it positive by replacing both a and d' by their negatives. This 

replacement will not change Using the circuit shown in Eq. ()101a|) . it is easy to 

(2) 

prove that is odd in both d and a'. 
Define 

= I [fee) I , ri = \ [abbc) \ = \ [ab)b ■ c\ . (131) 

To begin, we will assume that 7^ and rj 0. Later on, before ending the proof, 
we will remove these two constraints. 
If we define 

Xo = {a ■ c)[abb) ■ c , (132) 
Yo = {a-b){b-c)V , (133) 
(^;).=i,2,3=(c,[6'c'),&0, (134) 

a-W,3 = (a,^^,^), (135) 

i^j)j=i,2,3 = (a ■ bV'si3^, V'r], (d ■b)(b- c)(a ■ c )j , (136a) 
(«.W3 = (^.^^.^). (136b) 

■5/32 V 

(/ij)i=i,2,3 = (^-a' ■ c't], -(d ■ b){d' ■ c!)sp^, [abb) ■ cV'j , (137a) 
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(fi,W,3 = (^.^.c). (137b) 

then 

3 3 

= -Xo -tYo + J2 "jd'j^J + ' J2 . (138) 

3=1 j=l 

Define an angle (3 by 

cos(/5) = a' ■ c' , sm{(3)=V. (139) 
Define angles Pi, P2 ^ [0, tt) by 

cos(/5i) = a - b , sin(/3i) = |d x S| , (140) 

and 

cos(/?2) = b- c , sm{(32) = |& x c| . (141) 
Hence, [a6)/s/3^ = ^[6c)/s/32. One finds 

V = s/3Ac/3^\, (142) 

[d66c) 
V 

— [abbcc) 
V 

and 



^2 = -^2 , (143) 
■^3 = -a2, (144) 



[abb) ■ c = ^sp^Sfs^ . (145) 

At this point, it is easy re-express various quantities in terms of the principal 
parameters. Eq. (jl32p for Xo, Eq. (jl33|) for Yo, Eq. (jl36aj) for the Uj, and Eq. (jl37aj) for 
the iJ,j, yield, respectively, Ea. dT^, Eq.fTTTj) , Eq.(pH|l , and Ea.lfT^. 

We can also re-express Eqs. ()136b|) and ()137b|) for the vj and uj in terms of the 
principal parameters. One finds 

(%)i=i,2,3 = {92, -^^293, 9i) = {92, 93, 9i) , (146) 

and 

(%)i=i,2,3 = (-^6^2, ^3, ^1) = {92, 93, 9i) ■ (147) 
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Hence, for j = 1, 2, 3, 

V. = u - = g^^.-^ . (148) 

When S/32 or 77 vanish, Eq. ()135j) fails to define two of the vectors (jj, Eq. ()136bj) 
fails to define on or two of the vectors Vj, and Eq. ()137b|) fails to define on or two of 
the vectors Uj. If S/32 = 0, the proof survives if we define (^j)j=i,2,3 to be any RHON 
basis such that gi = c and §2 -L span{a,b,c). Then define the Uj and vj vectors 
in accordance with Eq. p48|) . If 77 = but S/jj 7^ 0, define the uj and Vj vectors in 
accordance with Eq. ()148p . 
QED 

Suppose we are given a matrix which is known to be the LO-RHS invariant 

(2'] 

^3 of a quantum circuit with 2-qubits and 3 DC-NOTs. Furthermore, we are asked 
to extract from this matrix values (non-unique ones) for a, b, c, a', b' and c'. Next 
we will give an algorithm for accomplishing this task. We will call it our "Algorithm 
for Diagonalizing Q^^" . The algorithm first expresses Q^^ in term of its principal 
parameters. Then it solves for a, b, c, a', b' and c' in terms of these parameters. 

(2) 

Algorithm for Diagonalizing ^3 : 



1. Set A3.. = |Re[tr(e;f )] and A3. = |lm[tr(6;f )]. Set A = - tr(6;f ), A, 



3r 



(A + At)/2 and Asi = (A - At)/(2z). Hence, ^^'^ = A3. + ^3. + A3, + lA^i, 
where X^r, As. are real scalars, and A^r, A^i are traceless Hermitian matrices. 



2. Set Xo = —Xsr and Fq = —A 



3j- 



Do a simultaneous SVD of A3, and Ag^. This decomposition is possible since we 
have shown previously that [A3^,A3j2. and [A3^,A3J/j are both zero. The de- 
composition yields orthogonal matrices U, V and real diagonal matrices D^r, D^i 
such that 



For i = 1, 2, 3, set 



Set 



A^, = UDsrV , A^, = UDs^V' . (149) 



[g[,g'„g',] = U , [^1,^2,^3]=^- (151) 
4. Set ^ = sign{fi3h'2) ■ Set ^2 = Calculate P from 

cb = e , ^° , s/3 = e , . (152) 
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If |c/3| > Isjsl, set 



C/3iC/32 



1 

C/3 



-/i2 

-6/^1 



(153) 



On the other hand, if \sp\ > Icpl, set 





C/3iS/32 


1 


Yo 


^1 


■S/3iC/32 


-5/315/32 






^/^3 _ 



(154) 



5. At this point, we know the four quantities Cjj^cp^, sp^cp^, C(j^S(j^, and sp^Sf^^. 
Calculate /3i ± [32 from 

cos(/3i ± [32) = cp^cp^ =F Sfs^Sfs^ , (155a) 

and 

sin(/?i ± /?2) = Sf}^cp2 ± C/3iS/32 • (155b) 
Calculate {Pi, P2) from /3i ± ^2- 

6. At this point, S/3jS/32 is guaranteed to be positive, but there is not guarantee that 
S/3j and s^2 individually positive (they may both be negative). Furthermore, 
at this point there is no guarantee that ^2 = sign{cp.^). These disagreements 
with the assumptions of our parameterization can be fixed as follows. If sp^ < 0, 
replace [3i and /?2 by their negatives, and replace {g[, §2,1^1,1^2, f^i, f^2) each by 
its negative. If ^2C/32 < 0, replace Pi ^ n — Pi and P2 n — P2, and replace 
{g[, §2, z/i, z/2, fii, /i2) each by its negative. 

7. Calculate a, b, c, a', b', c' from: 



c = 9i ( & = g[ 

b = Cfsji + S/32^3 ^ \ b' = g'^ . (156) 

a = cos(/32 - ^2pi)9i + sin (^2 - 6A)^3 [ a' = cpg[ + S/3^2 

Note the ^2's in the expression for a. The reason for these ^2's is that in order to 
obey —^^2 = +1, one must define the sign of the angle Pi differently depending 
on whether cp^ is positive or negative. (See FiglHI) 



Theorem 23 For any j G {1,2,3}, 

= XoY, . (157) 
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= c 



8 2 out 



' 1 = C S ^out 



C.>0 



<0 



Figure 9: Sign of (3i is defined differently depending on whether cp^ is positive or 
negative. 



Ifi,j,k are 3 distinct element of {1,2,3} , then 
and 



:i58) 



ly-Uj = -Yofik ■ (159) 

proof: 

Follows from the definitions Ea. lfT^ for Xo, Ea. (fT77|l for Yo, Ea. dT^ for the 
i/j, and Eq. (|129p for the /ij. 
QED 

Define 11 to be the permutation matrix that corresponds to the permutation 
map 7r() used above. Thus, 



n 



1 

1 

1 



(160) 



If (sj)j=i,2,3 denotes the standard basis, define matrices and by 



(161a) 
(161b) 



and 



M„ 



diag{vi,V2,^3)^ ■ 



(162a) 
(162b) 
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Note that A^^ ( given by Eq. p24p ) becomes My and (given by Eq. ()125p ) becomes 
when the bases (5'j)j=i,2,3 and (5'j)j=i,2,3 are both rotated into the standard basis. 



Theorem 24 



M^Mt = MM = XoYo , 



and 



proof: 
QED 



(Mj)2 = tr(M,) - M, . 



Follows from Theorem l23l 



8.4 Invariant for Circuits with 4 DC-NOTs 

[ ckt_invar4.m ] 

This part of our program is dedicated to the letters 
Theorem 25 




where 



where 



j 

Ki = Yoa^,j - a-, J: - a^,^^ + AY , 

^=^N ia'j ■ d') [g^(j)d) , y = (x)^^^ , 
j 

x' = ^f^j{9n{j) ■ d)[9jd') , y' = {x')f,^u , 



(163) 
(164) 

(165) 



(166a) 
(166b) 

(167) 

(168) 
(169) 

(170) 
(171) 
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AF = (AX) 



(172) 



where any variables not already defined in the statement of this theorem are defined 
in Section{E~M 

proof: 



-Qy- 










-(d)- 











An explicit expression for ^3 was given in SectionlHT 
the effect of DC-NOT similarity transformations. 
QED 



(173a) 
(173b) 

Eq. (P7|) shows how to calculate 



Theorem 26 When the bases (^j)j=i,2,3 O'lT'd (^j)j=i,2,3 o'^e both taken to be the stan- 
dard basis, then the quantities X^r, x, y, x' , y' , AX and AY (all defined in 
Theorem \25]} can be expressed in terms of the matrices M^^ and the vectors d,d' 
as follows: 



X^r — —d '^M^d , X^i — (A4r)i/- 
^ = [M^d', d) , y= (x)^^^ , 



x' = [M^d,d') , y'={x') 



(174) 
(175) 
(176) 



AX = d'(F{d '^mJ) - mJ(F - d'd '^M^ + , AY = (AX)^^^ . (177) 
proof: 

Just algebra. 

QED 

Theorem 27 See Fig\T^ 

Mjy' = YoX, M^x = Xoy', (178a) 
Mlx' = X,y, M,y = Yj'. (178b) 
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proof: 

Just algebra. 

QED 




Figure 10: Various vectors and what they are mapped into (up to a scalar factor) by 
and My. Since M'^ and Mj^M^ are both proportional to the identity matrix, 
one can replace by Mj and My by Mj in this figure if one also reverses the 
direction of the mapping arrows. 



9 Identities for Circuits with 2 Qubits 

This section deals with 2-qubit circuits, whereas Section [TUl deals with 3-qubit ones. 
In this section, with its numerous subsections, we start to reap the benefits of all our 
preceding hard work. The combination of dressed CNOTs and the LO-RHS invariant 
proves to be very useful. We find simple-to-check necessary and sufficient conditions 
for the reduction of a quantum circuit with j CNOTs to fewer CNOTs, where j = 2, 3. 
Plus we show how to express circuits with 1 or 2 controUed-U's as circuits with 2 or 
fewer CNOTs. Plus we show how to open and close a breach, a procedure that can 
reduce any 4-CNOT circuit to a 3-CNOT one. 

9.1 Reducing 2 DC-NOTs 

9.1.1 2 to 2 DC-NOTs (Angle Swapping) 

[ swap_aiigles .m, test_swang .m ] 

In this section we consider a circuit with 2 DC-NOTs acting on 2 qubits, and show 

(2) 

that a symmetry in Q2 allows one to swap certain angles without changing the effect 
of the circuit (up to LO-RHS). 
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As motivation for the main theorem of this section (the Angle Swapping The- 
orem), we present the next theorem. The next theorem shows that the target and 
control qubits of a controUed-U can be exchanged. 

Theorem 28 For any eR, 







g «2'^a 






-4>- 





(179) 



proof: 



g-4^a(0)e+^f'"S'(l) 



Jdaa (0) 1 ny (1) g-i|aa (0) g+i (1) 



(180a) 

(180b) 
(180c) 



QED 

The previous theorem immediately implies the next one, which states that we 
can "swap a breach" between two qubits. 





(181) 



Theorem 29 (Swapping a breach) Suppose 

For any C, it is possible to find an 71 such that C TZ. 
proof: 

Define 6 to be the angle between p' and q', and b the direction of p' x q'. Then 
p' ■ q' = cos{9) and p' x q' = sm{9)b' so ap'ag' = e^^'^y . Thus, 




(182) 



Given a unit vector a and an angle ^, we can always find (non-unique) unit 
vectors p and q such that angle{p, q) = 9, and pxq points along d. Then p - g = cos(6') 
and p X q = sin (6') a so apag — e*^°'". It follows that 




(183) 
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Now apply Theorem OHl to Eqs. ()182|) and 
QED 

Is it possible to swap a foil instead of a breach? Yes it is. In fact, one can 
swap any angle, as the following theorem shows. 



Theorem 30 (Angle Swapping) Let 
C = 




n 




:i84) 



For any L, it is possible to find an TZ such that C ~/j TZ and such that angle{b, a) 
angle{b'j,d'j) and angle(b',a') = angle{bf,df). 



proof: 



As proven in 

Section E21 'C^^^ 

can be parameterized as follows: 
















h 








(185) 



where a, a G M and where (/j)j=i,2,3 and (/j)j=i,2,3 are two RHON bases such that 



6 = /l , a = Cafl - Sah 



(186) 



and 



y = f[. a' = c^,f[-s^,f\. (187) 

can be parameterized in the same way as C^'^\ but with the replacements 
a af, a' -> a'j, fj ^ (/,)/, and f- -> {f))f. 

Our goal is to construct an IZ such that C IZ. Such an IZ, if it exists, must 
satisfy = We will use the positive sign. In light of Eq. (f7^ . this gives 

From the symmetrical form of the parameterized expressions for and 1Z^'^\ 
it is clear that these two invariants are equal if their principal parameters are related 
in the following way: 



/s/ — fl , flf — /s , f2f — — /2 



(189) 
(190) 
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and 



/s/ ~ fi ' fif ~ /s 5 /2/ ~ /2 • (191) 
These relations between the principal parameters of C^'^^ and 7^'-^^ imply that 




(192a) 



(192b) 



and 



7^ 




(193a) 



(193b) 



(Eq. ()193b|l is valid only if Sa and s^/ are both non-zero, whereas Eq. p93a|l is always 
valid. Theorem |2ni corresponds to the case = 0.) 

We are done proving the theorem, but we will go one step further, and give 
the value of the local operations U', U G SU{2) such that 

£ = 7^(t/'t f/t) . (194) 

When fi = fi = X and fs = fs = the right-hand sides of Eqs. ()192a|) and p93ap 
appear in Theorem ^2 It follows from Theorem ^2 and Eq. (^U)) that 



U 



(195) 



QED 



9.1.2 2 to 1 DC-NOTs 

In this section, we give necessary and sufficient conditions for a circuit with 2 DC- 
NOTs acting on 2 qubits to reduce to 1 DC-NOT. 
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o // o 

O -L O 



o ± cp 

o // 6 



Figure 11: All circuits with 2 DC-NOTs that reduce to 1 DC-NOT. 
Theorem 31 Suppose 



C 




n 




(196) 



d and 



For any C, it is possible to find an TZ such that L ~r IZ if and only if 
b' ± d') or (b J- d and V || d' ). See Fia\ll\ 

proof: 



Suppose b ^- d and b' \\ d' (the other case is analogous). When 6 ± a, 

^s(0)"'"'^^Va(0)"^'(^) = K,,(0)]'^^'(i) . (197) 
The last equation can be expressed diagrammatically as 







(b xdj 






1 


J 


L 



(198) 



Thus, when 6 ± d and b' = d' , C reduces to a single DC-NOT. More generally, 
d' = ±b'. Let Cnew be a new circuit obtained by replacing in C: d' by its negative if 
d' = —b'. By virtue of Eq.(j2SI), C = Cnewih ® U), where U G U{2). If Cnew ~/? T^new, 
then C ~R TZnew 



£ 7^ so £(2) = ±7^(2). In light of Eq.jZH), this gives 



It follows that 



where 



X2r + Asr + iA2i = ±i(Ta'^,aj 



(199) 



(200) 



= (d • b){a ■ b') 



(201) 
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A2r — ^ [a'b'b'), [abb) ' 



(202) 



A2. = a ■ ba^,^y f^ + a' ■ .^g . (203) 

= SO a ■ 6 = or a' • 6' = 0. Assume the former (the other case is analogous). 
Then a _L 6. = and a ■ 6 = so [d'b'b') = 0, which in turn implies that d' || b'. 
QED 



9.1.3 2 to DC-NOTs 

In this section, we give necessary and sufficient conditions for a circuit with 2 DC- 
NOTs acting on 2 qubits to reduce to zero DC-NOTs (i.e., to merely local operations). 

Q // O 



6 // 6 



Figure 12: All circuits with 2 DC-NOTs that reduce to DC-NOTs. 



Theorem 32 Suppose 



C 




(204) 



For any C, C 1 if and only if d \\ b and d' \\ b' . See FiaH^ 
proof: 

When d = b and d' = b', C equals 1. More generally, d = ±6 and d' = ±b'. 
Let Cnew be a new circuit obtained by replacing in £: (l)d by its negative if d = —b, 
(2)d' by its negative if d' = —b'. By virtue of Eq. ip^ . C = CnewiU' ® U), where 
f/', U E U{2). If Cnew 1, then C 1. 



C^rIso £(2) = ±1. In light of Eq.dni), this gives 

^^£(2) = ±1 . 

It follows that 

Asr + Asr + iK2i = ±1 • 

Thus = [d ■ b){d' ■ b') = ±1, which implies d || 6 and d' || b' . 
QED 



(205) 



(206) 
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9.2 Reducing 3 DC-NOTs 
9.2.1 3 to 2 DC-NOTs 

[ dr_3to2.m, test_dr_3to2 .m ] 

In this section, we give necessary and sufficient conditions for a circuit with 3 DC- 
NOTs acting on 2 qubits to reduce to 2 DC-NOTs. 

The constraint [abb) ■ c = shows up below. The field of Spherical Geometry 
sheds some light on this constraint. If we connect the points d, b, c by mayor-circle 
arcs on the unit sphere, then we get what is called a spherical triangle, [abb) ■ c = 
if and only if this spherical triangle has a right angle at vertex 6. (See Fig^^for an 
example of [d'6'6') ■ c' = 0.) 



right angle at * 


c 


p c 















6 6 6 



right angle at * 

Figure 13: All circuits with 3 DC-NOTs that reduce to 2 DC-NOTs. 



Theorem 33 Suppose 




n 




(207) 



For any C, it is possible to find an TZ such that C ~r TZ if and only if either [abb)-c = 
or [a'b'i)') ■ c' = 0. See Fig\T^ 

proof: 

Before we start the proof in earnest, let us restate some pertinent formulas 
taken from previous sections. 

From Section IHI21 we know that 



7^(2) 



\2r + A2r + i^2i 





f7 


fl 














Cq,' Sq, 



From Section ISTjI we know that 



(208a) 
(208b) 

(209) 
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where 



and 



A3, = [a'h'h') ■ c' [ahb) ■ c 



A3, = -(a ■ • c)V' - (a ■ h')(b' ■ c')V 



-{a' ■b'){d-b)c'c 



Asr = < 



A 



3i 



+ (a ■b)ib- c)[d'b'c')c' + id' ■ b'){b' ■ c^Yd[dbc) 
+ {d' ■ b')V[b'&)c^ + (a ■ b)V'&[bcf 

^ -[d'b'b'c'c')[dbbccf , 



+ (a ■ b)[d'b'c'c')[bccf + (a' • b')[b'c'c')[dbccf 
+ [dbb) ■ c[d'b'b'c')c^ + [d'b'b') ■ c'c'[dbbcf 



Now we begin the proof in earnest. 

C^rTZ so £(2) = ±7^(2). In light of Eq.^, this gives 



^£(2) = ±,2^(2) . 



It follows that 



= Xsr = [d'b'b') ■ c [dbb) ■ c . 
Thus, either [d'b'b') ■ c' or [dbb) ■ c. 
(^) 

Assume [d'b'b') ■ c' = 0. (The other case, [dbb) ■ c = 0, is analogous) 



a 





^ [ab') 



^ AAA 

A, [ a b' b') 
k' = — 



ki=b 



'I' 



Figure 14: Vectors and angles associated with bit-1 space spanned by d 



46 



It is convenient at this point to define a RHON basis {kj)j=i^2,3 for the 3d real 
space spanned by a',b',c'. Let sy = \ [a'b')\. If sxi 7^ 0, let 

If Sy = 0, define (fcj)j=i,2,3 to be any RHON basis such that k[ = b' and k2 i 
perpendicular to span(b', c'). Let 0' = angle{&, k'^). Since [a'b'b') ■ c' = 0, 



IS 



a' = cyk[ — syk'2^ , y = k[ , c' = s^/k[ + c^p^k'^ . (217) 

Eas. fl^TT)!) and (I^TTI) are illustrated in FiglTl 

Our goal is to construct an 71 such that £ 7^. Such an 71 must satisfy 
C^"^^ = ±7^'-^^ We will use the positive sign. In light of Eq. (f7^ . this gives 

t'C^'^ = ^27^(2) . (218) 

It follows that: 

\2r = -Xsi , (219a) 

= A3. , (219b) 

A2. = -Asi , (219c) 

A2i = Asr . (219d) 
By evaluating Eq. (j219a|) . we get 

Ca'Ca = {a-b){b-c)V + {a' ■b'){b' ■c')V (220a) 
= {a ■ b){b ■ c)syc^/ + cys^'V . (220b) 

Eq. ()219b|) is satisfied since [d'b'b') ■ c' = by assumption. 
By evaluating Eq. ()219c|) . we get 

f -{a -^Id' b'c'c')[bccf 

-Sa'Saflfl = I -{d' ■b')[b'c'c')[dbccf ■ (221) 

[ -[dbb)-c[d'b'b'c')c^ 



Define h by 



+syStf,'{d ■ b)[bcc) 

—cyC(p'[dbcc) ■ (222) 

+sy[dbb) ■ cc^ 
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If sy + and \h\ ^ 0, let 



[a'h'h'c') 



f - ^ 



(223) 



If \h\ = 0, set Sa'Sa = and choose any unit vectors for /2 and If 7^ but 
Sx' = 0, keep Eg . (12231) for s^'Sa and /2 but use /2 = ^ c'. 
By evaluating Eq. ()219d|) . we get 



= c'vj + k'2V2 



where 



and 



vi = —cx'{a ■ b)c + cx'S(i,i[dbc) + sx'C(f,/{a ■ b)[bc) 



(224) 



(225a) 



V2 = sx'Sff,'{a ■ b){b ■ c)c — cyc^'Vc + sx\abbcc) . (225b) 

At this point, we can follow from step El to the end of the Algorithm for 
Diagonalizing Q^"* that was given in Section |H21 This will yield values for a/, a'j, 6/, 
and b'j. 
QED 

Compared with the previous Theorem, the next theorem imposes more con- 
straints on £, and obtains a more constrained TZ. 

Theorem 34 Suppose 




, 7^ 




(226) 



Let A' = angle{d', b') and 0' = angle{c', a' x b'). For any C, if 

[ab'b') -0=0, 



and 



C0/(a ■ b)[db) - syCx'S^'b 



■c = 



(227a) 



(227b) 



then it is possible to find an TZ such that C TZ and such that b'j = c' . (Hence, c' 
"persists", from initial circuit C to final circuit TZ, as the bottom defining vector of 
the leftmost DC-NOT for both circuits. ) 
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proof: 



The (<^=) part of the proof of the previous theorem still applies. 

Using the definitions of Vi and V2 given by Eqs. (j225|) . it is not hard to show 



that 



v[v2 = 



c<^/(a ■ b)[db) - sycys^'b 



■c = . 



(228) 



Since Vi and V2 are orthogonal, the singular values and singular vectors of 
can be obtained simply by inspection of Eq. (j224j) . If | "Oil 7^ and \v2\ 7^ 0, then one 
can immediately set 



/s ~ ^2 



2 5 



1^2! 



F2I 



and 



r -r' f - 

Jl — C , J3 — T-^-T 



^1 



(229) 



(230) 



If \vi\ = but |'i;2| 7^ 0, choose /s = A x /2. If 7^ but |-U2| = 0, choose 
/i = /2 X /a- If = and |?T2| = 0, choose /i and /a to be any vectors that make 
(/j)j=i,2,3 a RHON basis. 
QED 



9.2.2 3 to 1 DC-NOTs 

In this section, we give necessary and sufficient conditions for a circuit with 3 DC- 
NOTs acting on 2 qubits to reduce to 1 DC-NOT. 

Theorem 35 Suppose 




n 




(231) 



Let V = d X b ■ c, and V' = d' x b' ■ & . For any C, it is possible to find an TZ such that 



■ ~iJ 


TZ if and only if one or more of the 




: ib\ 


a) and {V | d') 


Tib 


: (c| 


b) and {& \\ b') 


T^2a 


■■ {c' 


1 b' II a') andV = 


T2b 


: (c| 


b II d) and V = 




; d 1 


. span{b, c) and d' _L span{b', c'] 
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Figure 15: All circuits with 3 DC-NOTs that reduce to 1 DC-NOTs. The 8 circuits 
^±p,±p ^±q,±q defined by Eq. (|237p . 



T4 



span{d, h) and c' ± span{a', h') 



[ahh) -6 = and [a'h'h') ■ c' = 

\a-b\ \a'-b'\ \b-c\ \b'-c'\ 

sign {jlZr-r}! = " sign ^ ^ 



(232) 



(a-b)(b-c) 



(a'-6'){&'-c') 



proof: 



Consider a circuit of type T^b iT\a case is analogous). When c = h and c' = 6', it 
is obvious that a Ti^ circuit reduces to a single DC-NOT. More generally, c = ±6 and 
c' = ±6'. Let Cnew be a new circuit obtained by replacing in £: (l)c by its negative 
if c = -6, (2)c' by its negative if c' = -6'. By virtue of Eq.(I2SI), £ = ® U)Cnew, 
where f/', f/ G f/(2). If 7^„e^, then C ~r (?7' ® f/)7^„e^(f/'t ® W). 

Now consider a circuit of type (T2b case is analogous). Note that when 

V = 0, 



(233) 



so 





(234) 
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Thus, when a' = h' = c', a circuit reduces to a single DC-NOT. More generally, 
a! = ±b' and c' = ±6'. Let Cnew be a new circuit obtained by replacing in C: 
(l)a' by its negative if a' = —b', (2)c' by its negative if c' = —b'. By virtue of 
Eq.(I2a), £ = (/2 ® U)Cnew{l2 ® V^), where [/, G ?7(2). If Cm^n 7^„e«„ then 

(/2®t/)7^„e«,(/2®^/^). 

Circuits of type (T^f, case is analogous) reduce to a single DC-NOT by 
virtue of Theorem ITUl 

Now consider a circuit of type T4. For any Wi,W2 G {x,y,z} and ^ G M, let 
pI,^ ^2 ^'^d ^2 be defined as in Eq. . Because of the first line of Eq. ()232|) , one 
can choose a special coordinate system for bit such that c pf^, b x, a q^y, 
and a special coordinate system for bit 1 such that c' pf^, b' ^ x, a' ^ q^'y. See 

FigHni c, b, a and c', b', a' are portrayed in FigHni when (A;j)j=i,2,3 and {k'j)j=i^2,3 are 
the standard basis. In the special coordinate systems, the first hne of Eq. ()232|) is 
satisfied by construction. The second line of Eq. ()232|) becomes 



tan A I 



tan A' I and |tan( 



tan( 



(235) 



and the third line 



sign 



tan A 



sign 



/tan A' 
\tan (f)' 



(236) 



In general, Ea. (0!?K|l is satisfied iff A' G {±A, vr ± A} + 27rZ and 0' G {±0, vr ± 
0} + 27rZ. This gives 16 sign possibilities, but only 8 of them satisfy Eq. ()236|) . Indeed, 
let C^p,±p and C^q^^q denote the following 8 circuits: 



[Pi 



fP,Q 

'-"(-l)'"r,(-l)"r 



\§xy) 



[X 



-IV 



(237) 



where r G {p, q} and m,n & Bool. The following 4x4 matrix has rows labeled by 
the 4 possible values of 0', and columns labeled by the 4 possible values of A'. As its 
(0', A') entry, the matrix has: the T4 circuit implied by that value of (0', A'), if such 
a circuit exists, or an x if none exists. 



0'=i,A' =^ 


A 


TT- A 


TT + A 


-A 




X 




X 




TT — 




X 


r>P,q 
^-q-q 


X 


vr + 


X 


r>P,q 


X 


nP,q 


-0 


rp,q 


X 


np,q 
^q-q 


X 



(238) 



In conclusion, the 3 lines of Eq. ()232|1 imply, in the special coordinate systems, a circuit 
of type Eq.jSSZI)- 
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For (ditto, for C^'^p), there exists an 71 such that C TZ by virtue of 
Eq.(jSni) (ditto, Eq.(jnni))- The other 6 circuits of table Eq. ()238|) can be handled as 
follows. Let Cnew bc a new circuit obtained by replacing in £: (l)A' by A' — vr if 
A' = vr ± A mod (27r), (2)0' by 0' - vr if 0' = vr ± mod (27r). By virtue of Ea. ^ . 
C={U' ® U)CnewiV' (g) V) where U', U, V, V G f/(2), and where Cnew is of type CP'I 
or CP'l If Cnew ~K nnew, then C ~R (f/' ® f/)7^„eu,(^7't ® f/t). 



£ 7^ so £(2) = ±7^(2). In light of Eq.jZH), this gives 



It follows that 



where 



Asr + ^Asi + A.> + iA 



(239) 



(240) 



Asr = [d'b'b') ■ c [dbb) ■ c 



(241) 



\3i 



-(a-6)(6-c)V'- (a'- &')(&' -cOV 



(242) 



and 



A 



3r 



-(a' -6') (a •&)£'£ 
+ (a ■b)(b- c)[d'b'c')c^ + id' ■ b'){b' ■ &)&[dbcf 
+ (a' ■ b')V[b'&)c^ + (a ■ 6)V'c'[6c)^ 



-[d'b'b' c'c')[dbbcc) 



T 



(243) 



A 



3i 



+ {d ■ b)[d'b'&&)[bccf + (a' • b')[b'c'&)[dbccf 
+ [dbb) ■ c[d'b'b'c')c^ + [d'b'b') ■ c'c'[dbbc)'^ 



(244) 



b' if c'). 



First assume that there are no breaches in C (i.e., d ^ b, b ^ c, d' ^b', 
Note that 



[dbi)) -0 = and [d'b'b') ■ c' = . (245) 

This is why. From As^ = 0, we must have either [dbb) ■ c = or [d'b'b') ■ c' = 0. But 
if one of these holds, then the other one follows. Indeed, assume [dbb) ■ c = 0. Since 

also [db) ^ 0, it follows that [dbbc) 0. From A^i = 0, c'^A^i = [d'b'b') ■ c'[dbbc)^ = 
so [d'b'b') ■ c' = 0. By an analogous argument, assuming [d'b'b') ■ c' = leads to 
[abb) ■ c = 0. 
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Next note that 



a ■ b = a' ■ b' = =^ a ■ c = d' ■ c' = 



(246a) 



and 



Eqs. 



c-b = c' -b' = c-d = c' -a = . (246b) 
become obvious if one uses the BAG minus CAB identity to expand Eqs. ( 




Figure 16: Vectors and angles associated with bit-0 space spanned by a, b, c. Vectors 
and angles associated with bit-1 space spanned by d',b',c'. 



It is convenient at this point to define a RHON basis {k'j)j=i^2,3 for the 3d real 
space spanned by d',b',c'. Let sx' = \ [d'b')\. If sy ^ 0, let 

(247, 

If sx' = 0, let (%)j=i,2,3 be any RHON basis such that k[ = b' and k'^ is perpendicular 
to span(b', c'). Let cp' = angle{c', k'^). Since [d'b'b') ■ c' = 0, 

d' = cyk'i — Syk2, , b' = k[ , c' = S(j,ik[ + Cif,>k'^ . (248) 

Eas. pmi and (1^ are illustrated in FigUHl 

Use the previous paragraph with all the primes removed to define angles A, 
and a RHON basis {kj)j=i^2,3 for the 3d real space spanned by a, b, c. 

When expressed in terms of A, A', and 0', the constraint X^i = reduces to 

-[sx'Cxc^'S^ + ca'Sas^/c^] = . (249) 
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Likewise, the constraint Asj = reduces to 

T 

-[sx'Cxs^ic^ + cx'Sxc^'S^][k'^c')[k2c) = . (250) 
Eqs. fl249|) and ()250|) imply the following system of 2 equations: 

This system of equations can also be rewritten in the form: 

CySx SyCx 
_ SyCx Cx'Sx 

For Eq. ()251|) . either (i)the solution is the zero vector, or (ii)the determinant 
of the coefficient matrix vanishes. (i)If the solution is zero, then syCx = cySx = 0. 
Since we are assuming no breaches, sy 7^ and sx 7^ 0, so we must have cx = cy = 0. 
By virtue of Eq. (j246ap . this means that the circuit must be of type T^a- (h) If the 
determinant is zero, then 

I tan0 I = I tan0' | . (253) 

We will pursue this possibility later on. 

Likewise, for Eq. ()252j) . either (i)the solution is the zero vector, or (ii)the de- 
terminant of the coefficient matrix vanishes. (i)If the solution is zero, then Sfp/c^ = 
c^'S^ = 0. Since we are assuming no breaches, 7^ and c^f,' 7^ 0, so we must have 
S(f,' = Stj) = 0. By virtue of Eq. ()246b|) . this means that the circuit must be of type T^i,. 
(ii) If the determinant is zero, then 



SyCx 
cySx 



. 



(251) 



(252) 



|tanA I = |tanA' I . (254) 

We will pursue this possibility later on. 

Suppose we assume that the circuit C is not of type T3. Then, we have shown 
that it must satisfy Eqs. ()253|) and ()254p . But these two equations are the second fine 
of Eq. ()232p . To prove that the circuit must be of type T4, it remains for us to prove 
that the third line of Eq. ()232|) also holds. This third line clearly follows from X^i = 0, 
where X^i is given by Eq. (j242j) . 

Next , assume that there is at least one breach in C Without loss of 
generality, assume there is a breach between d and b (i.e., d || 6). 

d II 6 implies that V = 0. 

The constraint As^ = reduces to 

{b ■ c)V = , (255) 
which implies that either 6 ■ c = or V = 0. The constraint A^i = reduces to 

[a'b'c'c')[bccf = , (256) 
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which imphes that either (i)6 || c or (ii)a' || h' or (iii)[a'6') || c'. (i)If h || c, then, by 
Ea. d^ .^ V = 0. (a II b II c) and V = so £ is of type Tsb. (ii)If a' \\ 6', then, since 
also a II 6, £ is of type Ti^. (iii)Suppose [a'b') \\ c'. Assume that d' [f 6' as the case 
when these two vectors are parallel has already been considered. It follows that the 
conditions T^b are satisfied. This is why. [d'b') \\ c' and a' ^b' imply that V ^ 0, and, 
therefore, by virtue of Eq. ()255|) . c _L 6. Now c _L 6 and a || 6 imply that c _L a. Thus, 
c _L span{b,d). Also, [a'b') \\ c' implies that c' _L span{b' ,d'). 
QED 



9.2.3 3 to DC-NOTs 

In this section, we give necessary and sufficient conditions for a circuit with 3 DC- 
NOTs acting on 2 qubits to reduce to zero DC-NOTs (i.e., to merely local operations). 



Q // O // O 



O -L O -L O 



O -L O -L 



O // O // O 



Figure 17: All circuits with 3 DC-NOTs that reduce to DC-NOTs. 



Theorem 36 Suppose 

C= III • (257) 




For any C, C 1 if and only if one of the following is true (see Fig\r/{) 
Ta : (c', b', a!) are mutually orthogonal, and (c || b || a) 
Tb : (c, 6, a) are mutually orthogonal, and (c' || b' || a!) 

proof: 

Consider a circuit of type Tb {Ta case is analogous). Note that when {c,b,a) 
are mutually orthogonal, 

o'cO'bO'a = ic ■ [ba) = ±i . (258) 
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Hence, 




— (±iY-' — 



(259) 



Thus, when a! = h' = c', a Tf, circuit reduces to zero DC-NOTs. More generally, 
a! = ±b' and c' = ±6'. Let Cnew be a new circuit obtained by replacing in C: (l)a' 
by its negative if a' = —b', (2)c' by its negative if c' = —b'. By virtue of Eq. ()25|l . 

£ = (/2 ® U)CneM2 ® V) whcrC U,V & U{2). If Cne^ 1, then C r^R 1. 



£ 1 so = ±1. In light of Eq.jTH), this gives 



It follows that 



where 



±i 



(260) 



(261) 



Asr = [ab'b') ■ c' [abb) ■ c 



(262) 



A3, = -(a ■b){b- c)V - {a ■ b')(b' ■ c')V 



(263) 



and 



A 



3r 



-(a' -6') (a ■&)£'£ 
+ (a • 6)(6 ■ c)[a'6'c')c^ + (a' ■ b'){b' ■ &)&[abc)' 
+ {a' ■ y)V[b'&)c^ + (a ■ b)V&[bcf 
-[d'b'b'c'c')[dbbccf , 



(264) 



A 



C '^AarC 



+ (a • 6) [a'6'c'c') [6cc)'^ + (a' • 6') [6'c'c') [abccf 
~ ^ +[a66) ■ c[d'b'b'c')c^ + [d'b'b') ■ c'c'ldbbcf 



(265) 



i.3rc = SO d' ■ 6' = or a • 6 = 0. Both can't be true at once or else we 
would have X^i = 0, which is false. Assume henceforth that d' ■ b' ^ and d ■ b = 
(the case d' ■ b' = and d ■ 6 7^ is analogous). When d ■ 6 = 0, |A3j| = 1 reduces to 
|(d' ■ b'){b' ■ c')V\ = 1, which immediately implies that (c' || b' || d'), and (c, 6, d) are 
mutually orthogonal. Thus, circuit C must be of type T^. 
QED 
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9.3 Reducing Cont rolled- [/'s 
9.3.1 One Controlled-f/ 

In this section, we show that any controlled-U can be expressed with just two CNOTs. 
This resuh was first proven by Barenco et aL in Ref. [H]. Their method of proof is 
long and opaque compared with the uhra simple proof given below. This attests to 
the benefits of using dressed CNOTs. 



Theorem 37 Let 9 G 



Suppose 



n 




(266) 



For any C, it is possible to find an TZ such that C = TZ. 
proof: 

Given a unit vector w and an angle 9, we can always find (non- unique) unit 
vectors h and a such that angle{b, a) = 9, and bxa points along w. Then b-a = cos{9) 
and b X a = sm{9)w so cxgaa = e*^*^*. 



[g*e..{o)]n(i) ^ [a-^{0)a,{0)]-W = ag(0)"(iVa(0)"(i) 



(267) 



QED 



9.3.2 Two Controlled-[/'s (The Deflation Identity) 

[ def late_dcnots .m, test_def late_dcnots .m ] 

In this section, we show that a product of two controlled-Us can be expressed with 
just two CNOTs. This "Deflation Identity" was first proven in Ref. JHl- Unlike the 
proof of Ref . |10j , the one below uses dressed CNOTs. 



Theorem 38 Let A e SU{2) and 9l,9r e M. Suppose 



A 



n 




(268) 



For any L, it is possible to find an TZ such that C ~r TZ. 
proof: 

Given a unit vector wl and an angle 9l, we can always find (non- unique) 
unit vectors d and c such that angle{d,c) = 9^, and d x c points along wl. Then 
d - c = cos{9l) and dx c = sm{9L)wL so cr^ac = c^^^^'^l . Likewise, given a unit vector 
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A A 

[td) 

A A A A 

[ attd ) 




A A A A 

^ _ [ atdd) 

3 = ?j 



h^ = d 



AAA 

A [ at d ) 

h. = 

1 ri 




h^ = d 



h' = 



A^ A^ 

[ad) 



h\ = 



AAA 

[ a d' d' ) 



Figure 18: Variables used in Theorem 



Wi^ and an angle we can always find (non- unique) unit vectors h and a such that 
GiGa = e^^^^^^R. We are free to rotate the vectors d and c (ditto, h and d) within 
the plane they initially span, as long as we don't change the angle between them. In 
particular, we can choose both c and h to lie along the line of intersection between the 
planes span{d, c) and span{b, d). In other words, we can always choose c = b. Call t 
their common value . It is now clear that, without loss of generality, we can replace 
Chy 

— Q) — © — © — ^ — 

X II X X II X ■ '^^^^'^ 
—® — ® — ® — 

Our goal is to construct an TZ such that C ^rTZ. Such an 7^, if it exists, must 
satisfy = ±.1Z^'^\ We will use the positive sign. In hght of Eq.(f71|). this gives 

= ^27^(2) . (270) 
Using the same calculational techniques that were used in Section |Sl one finds 



£(2) 



(a ■ ■ d) - (a' ■ d')[dit) ■ d 

[a ■ t)d'[idf - (a' ■ d')d'[diidf + [d' d' d')[diddf 
+ [dt) ■ d[d'd')(F 



(271) 
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From Section IHI21 we know that IZ^'^^ can be expressed as 



7^(2) 



2i 



We must have 
















n 







(a'- 


d') [dii) 


■ d , 



(272a) 
(272b) 



and 



Define 



A2r = -Sa'SafUl = -[«^) " d[d' d')d 



= I [a (i') I , ?7 = I [did) \ . 
If 7^ 0, Ea. (l?fill is satisfied by 

.... 7 \d'd'' 

Sa'Sa = [at) ■ d S^> , /2 



f2 = d 



(273) 



(274) 



(275) 



(276) 



If s^' = 0, choose SaSa' and /2 the same way, but choose /a to be any vector perpen- 
dicular to d'. 

If Sif,' 7^ and r] 0, define the following two RHON bases (illustrated in 
FigdD: 



Ch' 



[d'd'd') [d'd') 



and 



1,2,3 



[didi) - [diddi) , 



?7 1] 



(277) 



(278) 



If s<)<,/ = 0, pick (/ij)j=i,2,3 to be any RHON basis such that /ig = d! . If 77 = 0, pick 
(/ij)j=i,2.3 to be any RHON basis such that h2 = d. Define the following two angles 
(illustrated in FiglTHj): 



02 = angle{[diid) , h^) , (pi = angle{[id), h^) 



(279) 
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We must have 



A. 



2i 



(a ■ t)d'[td) + {a ■ d')d'[attd) — [a d' d')[dtdd) 



hi 



"■3 



—d ■ ts^^ + c^'S(j,2 —d ■ tC(j,^ + C0'C(^2 

-S^'T] 



(280a) 
(280b) 



At this point, we can follow from step E] to the end of the Algorithm for 
Diagonalizing Q^^ that was given in Section IH^ This will yield values for 

and b'jr. 
QED 



9.4 Opening and Closing a Breach 

[ breach. m, test_breach .m ] 

Once more unto the breach, dear friends, once more; Or close the wall up 
with our English dead! (from "King Henry V" by W. Shakespeare) 

In this section, we show how to "open and close a breach" in 2-qubit circuits. 
This is a procedure whereby one can reduce any 2-qubit circuit with 4 CNOTs into 
a circuit with 3 CNOTs. Applying this procedure repeatedly, one can reduce any 
2-qubit circuit with more than 3 CNOTs into a circuit with only 3 CNOTs. The 
fact that all 2-qubit circuits can be expressed with 3 (or fewer) CNOTs was first 
proven in Ref. [B]. Unlike the proof below, their proof was based on Cartan's KAK 
decomposition [7J. 

Theorem 39 (Opening a Breach) Suppose 



(281) 



n= ^ X II X ■ ^^^^^ 

— © — © — (^iy- 

For any C, it is possible to find an IZ such that C IZ. 
proof: 

We begin by inserting a "unit wedge" into C: 
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In Eq. (l283|l . i and f are auxiliary parameters whose values are still to be defined. 

Consider separately each half of the circuit in Eq. ()28Hj) . Our goal is to re- 
express each half as follows: 



and 




From Theoreml^ we know that Eq. ()284|) will be achieved if we constrain our auxiliary 
parameters by: 

[PrQrQr) -^' = 0, (286a) 

and 

[c,t>'^iPR ■ Qr)Pr X ^i? - ■'^yu^x'a^^'n^R] - ^=0 ■ (286b) 

Likewise, Eq. (j285|) will be achieved if we constrain our auxiliary parameters by the 
same pair of equations as Eqs. (j286p . but with R subscripts replaced by L subscripts. 
These 4 constraint equations can be used to solve for the 4 degrees of freedom con- 
tained in the auxiliary parameters t and t'. 
QED 

By a "unit wedge" we mean a circuit element which equals one. An analogous 
concept is a "partition of unity". If it equals one, why use it? Because it depends 
on new, auxiliary parameters, and, by merging the unit wedge with its surroundings, 
we get a new expression which contains the auxiliary parameters, but is functionally 
independent of them. We can then choose convenient values for the auxiliary param- 
eters. The net result is that we can transform the original circuit to a new one that 
performs exactly as the old one but appears different. 

Note that in Eq. (j283|) we used a unit wedge consisting of a single DC-NOT 
times itself. There was no a priori obvious reason why this unit wedge would lead us 
to a proof of the theorem. We could have chosen a unit wedge that provided more 
auxiliary parameters. For instance, we could have chosen a product of 3 DC-NOTs 
(times the inverse of the product). After all, 1 DC-NOT can express only a limited 
subset of all possible 2-qubit transformations whereas 3 DC-NOTs can be used to 
express any of them. For proving the above theorem, using a unit wedge with only 
1 DC-NOT turned out to be sufficient. But one can envisage this theorem proving 
technique being used elsewhere with more complicated unit wedges. 
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Suppose one starts with a circuit which, hke £ in Eg . (12811) . possesses 4 DC- 
NOTs. By the last theorem, one can "open a breach" in it; that is, transform it into 
a circuit which, hke TZ in Eq. ()282|) . possesses two adjacent oval nodes both carrying 
a t' . Then one can combine the two adjacent DC-NOTs with a t' node and obtain 
a controUed-U. Finally, one can use the Deflation Identity of Sec ]9.3.2l to express the 
just created controUed-U and an adjacent DC-NOT as a circuit with two CNOTs. 
The net effect of this procedure is to reduce any 2-qubit circuit with 4 CNOTs into 
one with 3 CNOTs. 



10 Identities for Circuits with 3 Qubits 
10.1 Pass- Through Identities 

In the following 3 subsections, we consider the following 3 "identities" (one subsection 
per identity): 





-O 



(287a) 



(287b) 



(287c) 



-O 



Note that in all 3 identities, the initial and final circuits both have the same number 
of DC-NOTs, acting on the same 3 qubits. In all 3 cases, we pass a DC-NOT (the 
mobile one) acting on qubits and 1 through another DC-NOT (the static one) acting 
on qubits and 2. Thus, the mobile and static DC-NOTs both act on qubit 0, but 
the second qubit on which they act differs. We will refer to Eq. ()287a|) . Eq. ()287b|) . 
and Eq. ()287c|) as the Pass-Through Identities 1,2, and 3, respectively. In the initial 
circuit of Pass-Through Identity ra, the mobile DC-NOT is part of a group of n 
adjacent DC-NOTs acting on qubits and 1. 

The Pass-Through Identities Eqs. fl287p do not, per se, change the number of 
DC-NOTs. However, in some situations, they can be used to reduce the number of 
DC-NOTs. For example. 



-o- 



-o- 



'R 



(288a) 
(288b) 



-O 
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In Eq. ()288a|) . there are initially 3 adjacent DC-NOTs on the LHS of the static DC- 
NOT. Using Pass-Through Identity 1, we produce 4 adjacent DC-NOTs on the LHS of 
the static DC-NOT. As shown in Section 1^31 these 4 adjacent DC-NOTs can always 
be reduced to 3 DC-NOTs. 



10.1.1 Pass-Through Identity 1 
Theorem 40 Suppose 



C 





It 


If 












II \ 




V 7 



n 



(289) 



For any C, it is possible to find an TZ such that C TZ if and only if a || h. 
proof: 

(<^=) Let dj = a! and h'^ = b". Clearly, if a = 6, then C = TZ. More generally, a = ±b. 
Let Cnew be a new circuit obtained by replacing in C: a by its negative if d = —b. 
By virtue of Eq.(j2SI), C = Cnewih ®U ® I2) where U G U{2). If Cnew T^new, then 

C TZnew 
(^) 

Using the same calculational techniques that were used in Section |Hl one finds 
^^'^ = « ■ K",aM + ^^i,a',[aS) , (290) 



and 



(291) 



C TZ implies that C^'^^ is proportional to TZ^'^K Therefore, a^^, j-^^ must vanish. 

Hence, [ab) = 0, which is implies a \\ b. 
QED 

10.1.2 Pass-Through Identity 2 
Theorem 41 Suppose 




TZ 



(292) 
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For any C, if there exists v such that 





then it is possible to find an TZ such that L TZ. 
proof: 

One has 



(293) 




(294a) 



(294b) 



(294c) 



In (a), we introduced a unit wedge. To go from (a) to (b), we passed half of that unit 
wedge across the "static" DC-NOT. Finally, to go from (b) to (c), we used Eg . (12931) . 
QED 

Note that Section ID. 2. 21 gives necessary and sufficient conditions for a 2-qubit 
circuit with 3 DC-NOTs to reduce to an equivalent circuit with 1 DC-NOT. Using 
those necessary and sufficient conditions, it is easy to check in any particular instance 
whether there exists a t' such that Eq. ()293p is satisfied. 

10.1.3 Pass-Through Identity 3 

[ passS.m, test_pass3.m ] 

Theorem 42 Suppose 



C 




, n 




(295) 



64 



For any C, it is possible to find an IZ such that C IZ. 
proof: 




(296a) 



(296b) 



(296c) 



In (a), we introduced a unit wedge. To go from (a) to (b), we passed half of that unit 
wedge across the "static" DC-NOT. Finally, to go from (b) to (c), we used Theorem 



QED 

The next theorem is used in the proof of Theorem 021 
Theorem 43 Suppose 




n 




(297) 



For any C{-), there exists a d' and an TZ such that C ~ij TZ. 
proof: 

Our goal is to find a d' and to construct an TZ such that C TZ. Such an 
TZ must satisfy C^'^^ = ±TZ^'^\ We will use the positive sign. In light of Eq. lfTH) . the 
following must be true: 



From Section IH31 we know that 



(298) 
(299) 
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where 



= -d '^M^d 



A4» = -d'^M^d, 



Air = Xod'(F + x'(F + d'x^ + AX 
An = Yod'cF - y'(F - dY + AY . 



(300) 
(301) 
(302) 
(303) 



The precise definitions of {Xo,Yo), {x,x' ,y,y'), (AX, AF), and {Mfj_,M^) in terms of 
[a, a'), {b,b'), {c,c'), and {d,d') are given in Section lOl 
From Section IHI21 we know that 



Tl^^^ = X2r+A2r + iA2i 



el rT 





f7 


fl 









K 








(304a) 
(304b) 



We must have 



and 



A2r = — , 

= , 

A2r = —A^r , 

A2i = —A^i . 



(305a) 
(305b) 
(305c) 

(305d) 



To begin, we will assume that Xo 7^ 0. Later on, before ending the proof, we 
will remove this assumption. 

By evaluating Eq. ()305ap . we get 



Ca'Ca = d ^Mud . 
By evaluating Eq. ()305b|) . we get 



= d ' M^d . 



(306) 



(307) 



Let d' be any unit vector that satisfies this equation. 
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By evaluating Eq. ()305cp . we get 



Sa'SafUl = -{Xod'd'' + x'd'' + d' x'' + AX) . 



(30^ 



For Eg . (131)81) to be true, the RHS of that equation must factor into the product of a 
column vector times a row vector: 



Let 



where 



_ el cT 



-Xo d! + 



x' 



d + 



. d' + f 



V2 



, /2 



X 



V2 



V'2 



d' + 



x' 



1 + 



iX' 



'\2 



(Xo) 



, m = iv2) 



omit primes 



(309) 



(310) 



(311) 



Note that since Eqs. ()308j) and ()309|) are both true, the following must be true: 



x'x^ 
~Y7 



AX . 



(312) 



Eq. |312|) can also be proven by expressing it in terms of (a, d'), {b,b'), {c,c'), and 

{d,d'). 

By evaluating Eq. ()305d|) . we get 



A2, = -{Yod'd' - y'd' - d'y' + AY) . 



(313) 



At this point, we can follow from step El to the end of the Algorithm for 
Diagonahzing Q^^ that was given in Section IH^ This will yield values for d/, dj, bf, 

and b'j-. 

Now assume Xq = 0. By Eq. ()312|) . either x' = or x = 0. When x' = and 
x 7^ (the case x' 7^ and x = is analogous), Eq . (13091) becomes 



-Sa'Saf 2/2 



X 



d' + j^]x\ 



(314) 



where ^ is defined as the obvious limit. Thus, we can set 



/ I ^1 ff 
V2\-^\ ) J2 



d' + — 

V2 



/2 = 4. 

Ixl 



(315) 
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li x' = X = 0, then Eq. ()309|) becomes —Sa'Sa fliH = O5 so we can set Sa'Sa = and 
define /2 and to be arbitrary unit vectors. 

Additional observations: 

Note that /2 "^A2j = imphes 

X' ■ y' = XoYo , (316a) 

and 

AY^x' = Xoy. (316b) 
Likewise, note that A2j/2 = imphes 

x-y = XoYo, (317a) 

and 

AYx = Xoy'. (317b) 

Eqs. fl316|) and ()317p can also be proven by expressing them in terms of (d, d'), (b, b'), 
(c, c'), and (d, d'). 

If \x\ and \x'\ are both non-zero, it is possible to introduce 2 RHON bases 
(/ij)j=i,2,3 and (/ij)j=i,2,3, defined as follows. Define and /i2 by 



/^2 = /2, /^2 = /2. (318) 



Define h'^ and /i3 by 



'/ x'Xq 



d' - 



^3 — , y — (^3)omit primes ; (319) 

^3 



/\2 



where 



x'X 

d — 



[X 



i\2 



1 + ^^"^ = , ?73 = (%) omit primes • (320) 



Define h\ and /ii by 



[x'd') sign(X,) 



omit primes ; 



(321) 



where 



V'l = \ = \x'\ , Vi = Wi) 



omit primes 



(322) 
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After some algebra, one can show that Eq. ()313p becomes 









K 

\x'\ 


y^[xd) sign(Xo) 
-[x'd'f /\Y[xd) 


-Yo 

y"^[x'd') sign(Xo) 



(323) 



The entries of the previous table can be expressed solely in terms of {d, d') and 
{M^,Miy). After some algebra, one finds that 



f[xd) = iMjd')-[M^d'J) 



y'^[x'd') = {M,d) ■ [M^d, d') , 



and 



QED 



x'd!fAY[xd) = d^MlM^Mld' 



(324) 
(325) 

(326) 
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